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ABSTRACT 

We aim at the construction of luminosity and mass functions for Galactic open clusters, based on integrated magnitudes and tidal 
masses. We also aim at studying the evolution of these functions, with the ultimate purpose of deriving the initial luminosity and mass 
distributions of star clusters, independent of model assumptions regarding the cluster mass-to-light ratio. Finally we aim at a new 
determination of the percentage of field stars that have originated in open clusters. 

The integrated magnitudes are computed from individual photometry of cluster members selected from the ASCC-2.5 catalogue. 
The cluster masses we assumed to be the estimated tidal mass recently published by us elsewhere. Analysis of the cluster brightness 
distribution as a function of apparent integrated magnitudes shows that the cluster sample drawn from the ASCC-2.5 is complete 
down to apparent integrated magnitude I v = 8, with 440 clusters and compact associations above this completeness limit. This, on 
average, corresponds to a completeness area in the solar neighbourhood with an effective radius of about 1 kpc. 
The observed luminosity function can be constructed in a range of absolute integrated magnitudes l Uv = [-10, -0.5] mag, i.e. about 

5 magnitudes deeper than in the most nearby galaxies. It increases linearly from the brightest limit to a turnover at about I Mv « -2.5. 
The slope of this linear portion is a = 0.41 ± 0.01, which agrees perfectly with the slope deduced for star cluster observations in 
nearby galaxies. The masses of the Galactic clusters span a range from a few M Q to logM c /M G » 5.5. The mass function of these 
clusters can be fit as a linear function with log mass for log M C /M Q > 2.5, and shows a broad maximum between log M C /M Q = 1.5 
and 2.5. For log M C /M Q > 2.5, the linear part of the upper cluster mass function has a slope a = 2.03 ± 0.05, again in agreement 
with data on extragalactic clusters. We regard this agreement as indirect evidence that the tidal masses for Galactic clusters and the 
luminosity-based masses for extragalactic clusters are on the same scale. 

Considering the evolution of the cluster mass function now reveals a slight but significant steepening of the slope with increasing age 
from a = 1.66 ± 0.14 at log t < 6.9 to a = 2.13 ± 0.08 at log / < 8.5. This indicates that open clusters are formed with a flatter (initial) 
mass distribution than the overall observed cluster mass distribution averaged over all ages. Interestingly, the luminosity function of 
open clusters does not show the same systematic steepening with age as the mass function does. 

We find that the initial mass function of open clusters (CIMF) has a two-segment structure with the slopes a = 1 .66 ± 0. 14 in the range 
logM r /M = 3.37 . . .4.93 and a = 0.82 ± 0.14 in the range log M C /M Q = 1.7.. . 3.37. The average mass of open clusters at birth is 
4.5 • 10 3 M G , which should be compared to the average observed mass of about 700 M Q . The average cluster formation rate derived 
from the comparison of initial and observed mass functions is u = 0.4 kpc~ 2 Myr~'. Multiplying by the age of the Galactic disc (T 
= 13 Gyr) the predicted surface density of Galactic disc field stars originating from dissolved open clusters amounts to 22M Q pc~ 2 
which is about 40% of the total surface density of the Galactic disc in the solar neighbourhood. Thus, we conclude that almost half of 
all field stars were born in open clusters, a much higher fraction than previously thought. 

Key words. Galaxy: disk - Galaxy: open clusters and associations: general - solar neighbourhood - Galaxy: stellar content - 
Galaxies: star clusters 



1. Introduction 

Open clusters constitute an important part of a process trans- 
forming gas and dust into stars. They are observed as the most 
prominent parts in the regions of active star formation, or as trac- 
ers of the ceased star formation process in the general Galactic 
field. However, the role they are playing in this process has 
still not been fully understood. In spite of their prominence, 
there are indications that classical open clusters contribute 



only 10% or even less in put dWielenl Il97lb iMiller & Scald 
ll978UPiskunov et alj|2006l) to the total stellar population of the 
Galactic disc. This contradiction can be explained either by an 
early decay of a considerable fraction of newly formed star clus- 
ters ( see e.g. lTutukov|[T97l iKroupa etafj 120011: lLamers et ail 
120051) or by an insufficient knowledge of cluster formation statis- 
tics. In this context, one should note that the most important 
items of cluster formation like the distribution of cluster masses 
at birth (i.e., the initial mass function of star clusters) and the 
cluster formation rate are still poorly known. 
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Even in the close vicinity of the Sun, the only at- 
tempt to construct the luminosity function of open clusters 
( Ivan den Bergh & Lafontainel 19841) is based on a sample of 142 
clusters that, according to the authors, is to 2/3 complete within 
400 pc. A mass function of Galactic open clusters is not yet 
available. 

In contrast, extragalactic cluster populations have been ac- 
tively investigated, especially during the present decade. The lu- 
minosity functions of remote cluste rs are determined both in ac- 
tive and regular galaxies (see e.g. iLarsenl 120021; Ide Griis et alj 
120031: iGieles et all 120061 and references therein). For some 
galaxies, luminosity distributions of open clusters are con- 
verted to the mass functions via the as sumed theoretical mass- 
luminosity relation (Zhang & Fall 1999). 

Using extragalactic data for the study of cluster population 
gives certain advantages. When observing external galaxies, one 
is able to cover large areas of their surfaces and catch the bright- 
est objects, which are intrinsically rare. This is rather difficult in 
the case of our own Galaxy since they could be hidden by a dusty 
environment. Also, different selection effects can be taken into 
account better in an external galaxy because its clusters are lo- 
cated at about the same distance from the Sun. Moreover, studies 
of extragalactic clusters give a key to our understanding of the 
cluster formation in different environments. 

However, there are a number of problems related to studies 
based on extragalactic clusters. Samples of extragalactic clusters 
are restricted to relatively bright objects. There is a danger of a 
contamination of cluster samples by brightest stars and o f losing 
sparse objects in the strong background. For example, ILarsenl 
d2002l) finds that his sample in NGC6949 is contaminated to 
about 20% by bright stars. The major problem is the necessity 
of transforming the observed luminosities into cluster masses. 
The corresponding transformations are mainly based on the lu- 
minosity of a dozen of the most massive stars, whereas cluster 
masses are defined by numerous stars of lower masses that are 
below the observing limit. 

In our study of the local population of Galactic open clusters, 
we obtained data that now allow a reliable construction of their 
luminosity and mass functions. Our cluster sample contains 641 
open clusters and 9 associat ions identified in the all-sky com- 
piled catalogue ASCC-2.5 dKharchenkoll200 ll) . For each cluster 
we determined a homogeneous set of cluster parameters includ- 
ing, among others, its distance, age, and tidal mass. The sam- 
ple is found to be complete up to a distance of about 850 pc 
from the Sun (Pis kunov et al.l 120061) . This completeness limit 
corresponds to a dis tance modulus (V - My) of about 10-10.5 
dSchilbach et al.|[2006h and allow probes of the luminosity func- 
tion of open clusters that are much deeper than in other galaxies. 

In this paper we aim at constructing the luminosity and mass 
functions of the Galactic open clusters and at comparing them to 
the data in other galaxies in order to touch the issue of univer- 
sality of the cluster formation process in different environments. 
We also consider temporal variations of the two distributions. 
The ultimate purpose of this study is to construct the initial lu- 
minosity and mass distributions of star clusters independent of 
model assumptions on the M/L-ratio, to derive their parameters 
(the shape and limits), and to draw conclusions on the impact of 
these results on our understanding of the role of open clusters in 
the evolution of the stellar population of the Galactic disc. 

The paper has the following structure. In Sect. [2] we briefly 
describe the input data and give the definitions used throughout 
the paper. Sect.[3]is devoted to the construction of cluster the lu- 
minosity and mass functions. In Sect. [4] we examine how cluster 
mass and luminosity functions evolve with time and build the 



cluster initial mass function. In Sects. [5] and [6] we discuss and 
summarise the results. 



2. Data and definitions 

In this paper we use the results of the previous work based on 
the study of 650 nearby open clusters identified in the ASCC- 
2.5. The sample includes 520 kn own objects and 130 open clus- 
ters detected within our project (Kharchenko et al. 2005b). For 
each star projected on a cluster area, a membership probability 
was determined in an iterative process that takes spatial, pho- 
tometric, and proper motion distributions of stars into account 
withi n the corresponding area on the sky (Kharc henko et al.l 
2004). At the end of the iterations for each cluster, we ob- 
tained new coordinates of the cluster centre, the cluster size, the 
mean proper motion, the dista nce from the Sun, reddening, and 
age dKharchenko et al. 2005a b). These parameters were deter- 
mined with data on the most probable members, i.e., stars hav- 
ing both kinematic and photometric membership probabilities 
higher than 61%. The results are included i n the Catalogue of 
Open C luster Data (COCD) and its Extension dKharchenko et al.l 
\2005M . 

In iPiskunov et al.l d2007l) we published masses of 236 clus- 
ters of our sample estimated from tidal radii determined with a 
three-parameter fit of King's profiles to the observed density dis- 
tributions of cluster members. To obtain mass estimates for all 
clusters, we used these data to establish a relation between tidal 
radius and the observed semi-major axis of the apparent distri- 
bution of cluster members on the sky. The resulting tidal masses 
(i.e. calib rated tidal masses) f or 650 clusters have been pub- 
lished in Piskunov et al. (2008). Though calibrated tidal masses 
of individual clusters are less accurate, their distribution fits well 
the distribution of the tidal masses based on the direct fitting of 
the King model. Moreover, because the parameters are homo- 
geneous and more numerous, the calibrated masses are better 
suited to statistical investigations. In the following, we use the 
data on the calibrated tidal masses to study the luminosity and 
mass functions of the Galactic open clusters. Throughout the pa- 
per, masses are given in the units of solar masses. 

Our list of 650 Galactic open clusters identified in the ASCC- 
2.5 includes 9 cluster-like associations (e.g. Vel OB2, Sco OB4), 
which require further commentary. These are objects, which are 
classified both as clusters and associations in the literature (e.g. 
a Per = Per OB3). They are included in our sample because 
their observed properties (spatial structure, internal kinematics, 
and stellar content) make them indistinguishable from regular 
open clusters, although they have larger sizes and higher masses. 
Hereafter we refer to th em as cluster-like associations (see also 
iKharchenko et alj|2007l) . 

We consider two young clusters, NGC 869 (also known as 
h Per) and NGC 884 (x Per), as a single entity. They overlap 
in the projection on the sky, share a huge corona, and have a 
large number (more than 50%) of members in common, mak- 
ing an accurate determination of their individual tidal radii and 
masses rather difficult. Also, we exclude the cluster Mamajek 1 
(rj Chamaeleontis) for which we were able to identify only three 
members in the ASCC-2.5. According to Fig. [2] this cluster can 
be omitted without consequences for the results. Therefore, our 
final cluster sample includes 648 entities. 

As a measure of cluster brightness needed for the construc- 
tion of the luminosity function, we take the integrated magni- 
tudes of the clusters. Doing so, we can directly compare our 
results with findings based on observations of open clusters in 
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Fig. 1. Absolute integrated magnitude Im v profiles for 9 "tem- 
plate" clusters. AV is the magnitude difference between the 
brightest cluster star and fainter members. The vertical bars in- 
dicate individual stars. The contribution of stars fainter than 
AV = 7 to Im v is small enough to be neglected. As an illustra- 
tion of the reduction method applied to the observed integrated 
magnitudes of clusters, we show the profile of ASCC 51, which 
ends at AV a 5 mag (longer vertical bar). The extension (dashed 
curve) to AV = 7 mag is based on the template Mel 111, and it 
results in Aly = 0.007 mag. 

other galaxies. We define apparent and absolute integrated mag- 
nitudes Iy and Im v as 



/y = -2.5 log 
and 



A'l 



0.4 V, 



+ 10 



-0.4 M v 



Im v =I v -(V- M v ) , 

where N\ and V; are the number and the apparent magnitudes 
of the most probable cluster members, respectively, and (V - 
My) is the apparent distance modulus available in the COCD for 
each cluster. For "unseen" stars, Aly is a magnitude correction, 
introduced to make the computed Iy and Im v independent of the 
range of stellar magnitudes actually observed in a given cluster. 

To compute Aly, we used (Im v , AV) profiles of nine "tem- 
plate" clusters shown in Fig. Q] The templates were selected 
to satisfy two constraints: to be uniformly distributed over the 
brightest stellar magnitudes of the clusters and to represent a suf- 
ficiently long observed main sequence (MS). In Fig. Q] we show 
the profiles as a function of AV; = V, - VV (V; are individual 
magnitudes of the most probable members of a template and 
Vbr = min{V,}). At AV; = 0, the integrated magnitude Im v is 
identical to the absolute magnitude of the brightest cluster mem- 
ber. As one can see, stars that are seven or more magnitudes 
fainter than the brightest member (AV > 7) do not have a sig- 
nificant impact on Im v - On the other hand, the typical absolute 
magnitude of cluster stars at AV = 7 is brighter than My — 6 or 
m > 0.9 M in our sample. In this mass range, the mass func- 
tions of cluster stars are only weakly dependent on cluster age 
(Baumgardt & Makino 2003). This enables us to safely apply the 
above templates to the whole cluster sample. 



A correction Aly is computed for each cluster having a short 
MS as 



Ah = I f v n ' - 2.5 log 2 

v,=v f „ 



10 



-0.4 V, 



after normalization of the suitable template profile (see Fig. [TJ. 
Here V,,- and V/,„ are the brightest and faintest magnitudes of the 
most probable members in the cluster under consideration, V, are 
template magnitudes, and AV = 7 mag. The constant I v "' is com- 
puted assuming that the cluster and the matching template have 
equal integrated magnitudes at V/„ t . Since the integrated magni- 
tude is mainly defined by the brightest cluster members, Aly is 
relatively small: it is always less than 0.85 mag, with an average 
of 0.13 mag, and for 95% of the clusters Aly < 0.3. On the other 
hand, in 33 clusters the MS exceeds a length of 7 magnitudes. 
We truncated their profiles to AV = 7, for homogeneity. 

The other membership samples, i.e., samples including stars 
with a membership probability less than 61%, produce some- 
what brighter magnitudes (on average, by about half a magni- 
tude) but they are more strongly contaminated by field stars than 
is the sample of the most probable cluster members. 

Throughout the paper we use the following definitions of 
cluster luminosity and mass functions. If the AN is the num- 
ber of open clusters having absolute integrated magnitude in the 
range [Im v , Im v + &Im v ] and which are observed in an area AS 
in the Galactic disc, then the luminosity function </> is 



<I>(Im v ) = 



1 AN 



AS AI 



My 



i.e., the luminosity function is a surface density distribution of 
open clusters over the integrated magnitude. Similarly, the mass 
function rj is defined as a surface density distribution of open 
clusters over the logarithm of mass: 



V(M C ) 



1 



AN 



AS AlogM, 



Since the above distributions include clusters of all ages, it is rea- 
sonable to call 4>(Im v ) and r](M c ) the present-day luminosity and 
mass functions of clusters, or CPDLF and CPDMF, respectively. 
In the following, we also consider the luminosity and mass dis- 
tributions of clusters confined by some upper limit of their age 
t. We call them current luminosity/mass functions of clusters or 
simply cluster luminosity/mass functions with the abbreviations 
CLF and CMF To distinguish them from present-day distribu- 
tions we denote them as <p r (lM v ) and j],(M c ). Both functions are 
in fact cumulative with respect to age distributions. 

The aim of the present paper is the construction of an initial 
mass function of star clusters, which indeed describes the initial 
distribution, i.e. th e distribution after re-v irialisation after resid- 
ual gas expulsion (iKroupa & Boilv 2002). Hereafter, we denote 
the initial luminosity/mass distributions as CILF and CIMF. A 
formal definition of these functions is given in Sect|4] 

Finally, we must keep in mind the evolutionary status of open 
clusters included in our sample. Since cluster membership is 
based on the proper motion data mainly obtained in the opti- 
cal spectral range, we consider our sample as representative of 
optical clusters or "classical" open clusters. The embedded ob- 
jects (e.g., the nearby cluster NGC 1333) are not included in our 
statistics since their members would be fainter than the limiting 
magnitude of V ~ 11.5 of the ASCC-2.5 if observed in the op- 
tical. Therefore, we assume the beginning of the transparency 
phase after the removal of the bulk of the placental matter to be 
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Fig. 2. Distribution of 650 clusters of our sample in 
(log t , log M c )-diagram (dots). Open circles mark volume- 
limited sample of 257 clusters with a border of 850 pc. Crosses 
are 440 clusters from the magnitude-limited sample constructed 
in Sect. [3] Triangles mark two clusters omitted from the consid- 
eration (see Sect. |2j. The dashed line shows the upper age limit 
of the youngest subsample discussed in Sect|4] 

a starting point of the evolution of a classical open cluster. The 
corresponding age to is defined by the lowest age of our clusters, 
that is, about 4 Myr. 

3. Magnitude-limited sample of open clusters and 
the observed luminosity and mass functions 

As shown in Pisku nov et alj d2006l) . our sample can be consid- 
ered as volume-limited with a boundary of about 850 pc. This 
volume contains 2570 clusters, i.e. about one third of the total 
sample. However, the completeness distance varies if one con- 
siders clusters with different integrated magnitude. Very faint 
clusters can possibly be seen within slightly smaller distances. 
In contrast, very bright and massive clusters are seen at much 
greater distances but they are rare, and the limit of 850 pc is not 
sufficient for a realistic estimate of their average density. This 
presumption is supported by Fig.|2]where in an "age-mass" dia- 
gram we compare the distributions of all clusters with the clus- 
ters belonging to the volume-limited sample. For log t < 7.5, the 
volume-limited sample does not represent the actual mass dis- 
tribution of the young clusters and cannot be used for a realistic 
estimate of a mass function. The same is true for the distribution 
of integrated absolute magnitudes of young clusters. The prob- 
lem can be solved if we could extract a cluster sample that is 
limited by a certain integrated magnitude. In this case, we can 
determine the completeness areas for clusters of different intrin- 
sic brightnesses. 

In Fig. [3] we show the distribution of our clusters in bins of 
apparent integrated magnitude Iy. Starting from Iy w 3, the log- 
arithm of cluster numbers is linearly increasing with increasing 

1 This number does not include the two nearest open clusters UMa 
and the Hyades, for which we could not determine cluster membership 
via the standard pipeline and, therefore, were not able to obtain cluster 
masses and integrated magnitudes. 




0.0 i J. J. 1.1. J. i I i i i I i i i I i i i I i i i 
2 4 6 8 10 12 

I v , mag 

Fig. 3. Distribution of open clusters versus apparent integrated 
magnitude Iy computed from data on the most probable cluster 
members (crosses). The bars are corresponding Poisson errors, 
the solid line is a linear fit of the distribution and the vertical 
dotted line is the adopted completeness limit Iy. 

Iy and reaches a maximum at Iy = 8. We take Iy — 8 as the 
completeness magnitude of our sample. In total, 440 clusters are 
brighter than Iy. This increases the basis for the following statis- 
tical study by a factor of ss 1 .7 (compared to the previous number 
of 257). 

The cluster luminosity function (p is constructed as a function 
of integrated absolute magnitudes Im v j binned with a step size 
AIm v = 0.25, as the sum of partial surface densities 

1 1 

<f>(lM Y ,d = ~ry J] ~32 ■ 

v j j 

Here dj is the individual completeness distance corresponding to 
Iy = 8, n, is the number of clusters in the 2th magnitude bin Im v ,u 
and Ya n i - 440. Similarly, the mass function 77 is computed by 
using the logarithmic mass scale log M c j binned with a step size 
AlogM c = 0.15 

with v, the number of clusters in the mass bin log M c j and £ v, = 
440. The completeness distance dj for the jth cluster is given as 

log dj = 0.2 [l v - I MvJ - 3.1 Ej(B - V)] + 1 , 

where Ej(B - V) is the reddening for the cluster. 

In Fig. |4] we show the present-day luminosity function 
CPDLF constructed from 440 clusters. The CPDLF splits into 
a long quasi-linear portion at bright magnitudes, a turnover be- 
tween Im v = -3 and -2, and an apparent decrease at fainter 
magnitudes. As expected, for Galactic star clusters the CPDLF 
can be observed as much deeper than for clusters in other galax- 
ies. For the Large Magellanic Cloud (LMC), the faint limit is 
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Fig. 4. Luminosity function of Galactic open clusters based on 
440 local clusters brighter than the completeness limit Iy of the 
sample. The bars are Poisson errors, the dashed line shows a 
linear fit for the brighter part of the histogram (Im v < -2.5) 
where a is the corresponding slope. The arrow indicates the limit 
of integrated a bsolute magnitudes reached for open clusters in 
the LMC (see lLarseril2002t) . 

reached already at abo ut Im„ — -5 , and it is much brighter in 
more distant galaxies (Larsen 2002]). 

Since the luminosity function is obtained from the distri- 
bution of clusters within the respective completeness limit, the 
turnover can be considered as a real feature of the CPDLF in- 
dicating a decrease of cluster density at Im v > -2.5. The slope 
a on the left side of the turnover (Jm v < -2.5 ) is computed by 
fitting to the histogram (Fig.|4| a function 

log tp = b + a I Mv , 

where a = 0.41 + 0.01, b = 2.58 + 0.07. This corresponds to a 
slope of the cluster luminosity function a = 2.02 + 0.02 counted 
in luminosity intervals (for a compatibility with stellar luminos- 
ity function we adopt that dN/dLy °c L~ a ). The slope a is fre- 
quently used for the descripti on of the lu minosity function of ex- 
tragalactic clusters (see e.g. lLarsenll2002l) and can be expressed 
with the adopted definition as a — 2.5 a + 1. 

The only attempt to obtain the luminosity func- 
tion of Galactic open c luster s was undertaken by 
Ivan den B ergh & Lafont aind d!984l) . These authors con- 
structed the CPDL F from integrated magnitudes computed by 
Saga r et ail d 19831) for a sample of 142 clusters with accurate 
U BV photometry and proper motion membership. Their CPDLF 
was obtained for Im v < -2 and with a slope a = 0.2 mag -1 , 
which is flat compared to our result. The authors regard the 
sample to be 2/3 complete within 400 pc and strongly incom- 
plete outside this limit. However, no further details are given 
about the completeness issue (e.g. how it depends on cluster 
magnitude), and it is not possible to judge if their low value for 
the slope is a con sequence of the incompleteness. 

iLarsenl d2002l) studied the luminosity functions of open clus- 
ters in six nearby spiral galaxies constructed from HST archive 



Fig. 5. Mass function of Galactic open clusters based on 440 lo- 
cal clusters brighter than the completeness limit Iy of the sample. 
The bars are Poisson errors, the dashed line shows a linear fit for 
the mass rich part of the histogram (log M c > 2.5) where x is the 
corresponding slope. 

images, and in the LMC from literature data. Typically derived in 
the magnitude range [-10, —7], the CPDLF slope varies within 
a x 0.4 . . . 0.6 mag 1 . This result coincides with our finding for 
open clusters in the Milky Way. Further evidence for an agree- 
ment of the CPDLF slopes b etween Galactic an d extragalactic 
clusters can be found e.g., in lde Griis et al.l d2003l) . 

The other interesting feature in our CPDLF is a deficiency of 
clusters for Im v < -8 co mpared to the fitted power law. A similar 
effect was observed by ILarsenl d2002l) for extragalactic clusters. 

The corresponding present-day mass function CPDMF of 
440 clusters is shown in Fig. [5] We note similar details in the 
CPDMF that we have just observed for the CPDLF in Fig|5] One 
recognises a quasi-linear portion in rj(M c ) for logM c > 2.5, a 
turnover at about log M c = 2, and a decrease for smaller masses. 
Again, we stress that the turnover of the CPDMF can be con- 
sidered as a real feature, since the determination is based on 
unbiased cluster data. In the high-mass part of the histogram 
at log M c > 2.5, the CPDMF can be expressed in a canonical 
power-law form as 

log r] = log tj*-x log M c . (1) 

A fit of Eq. [T]to the observed mass function provides x — 1.01 + 
0.04 for the slope and log rf - 4.3 + 0.14 for the zero-point. 
A deficiency of clusters at the high-mass end (logM c ^ 4.0) 
is less pronounced than we observe for the luminosity function, 
and discrepancies to the fitted relation are almost within Poisson 
errors. 

Since all published data on mass functions of extragalac- 
tic clusters are based on the mass-luminosity relation used to 
convert the observed photometric or spectroscopic data into 
m asses, only a formal comparison with our result is possible. 
In Ide Griis et at] d2003l) . the slopes of the CPDMFs are listed 
for open clusters observed in four galaxies. Within a mass range 
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Fig. 6. Evolution of the mass function T],(logM c ) of Galactic 
open clusters. Different symbols mark samples with different up- 
per limits of cluster ages. The filled circles are for clusters with 
logf < 6.9, stars show CMF for log? < 7.9, and crosses indi- 
cate the CPDMF based on all 440 clusters (log f < 9.5). The 
bars are Poisson errors. The straight lines are the corresponding 
fits to linear parts of the mass functions at masses greater than 
logM c = 3.3 indicated by the vertical dotted line. The curve is 
the smoothed CPDMF at logM e < 3.3. The arrow indicates the 
lower mass limit reached for open clusters in the LMC. 



Fig. 7. Evolution of the slopes of the mass function j],(logM c ) 
at log M c > 3.3. The bars are Poisson errors. 



ity and mass functions. Whereas the cluster luminosity function 
evolves in coordination with the evolution of the most massive 
cluster members, the mass function evolves faster at its low-mass 
end where the dynamical time is short. 

Since the cluster luminosity function shows a more compli- 
cated evolution pattern, we start our consideration with the CMF. 



(10 3 , 10 6 M o ), a typical value of the CPDMF slope is a * 2. 
On the logarithmic scale, this corresponds to^ » 1 and agrees 
well with our results. We consider this result as indirect evidence 
of the coincidence of mass scales of Galactic and extragalactic 
clusters. 

4. Evolution of luminosity and mass functions; 
cluster IMF 

Since star clusters evolve by changing the basic parameters (like 
mass and integrated luminosity), their mass and luminosity func- 
tions also undergo evolutionary changes. For stars, the impact of 
their evolution on the stellar mass and l uminosity functions was 
revealed for the first time by ISalpeterl d 19551) . who developed 
a receipt for the construction of stellar initial luminosity/mass 
functions from present-day distributions. 

Compared to stars, the case of stellar clusters is rather com- 
plicated. The evolution of a cluster follows two independent time 
scales. The first one is the nuclear time scale that governs the 
evolution of stars. The nuclear scale is primarily responsible for 
changing the cluster's luminosity and is related to the luminos- 
ity function of clusters. The second time scale is defined by the 
dynamical evolution and, especially, by two-body-relaxation- 
driven mass loss from clusters. Therefore, it is responsible for 
the evolution of the cluster mass function. Moreover, the nuclear 
lifetime does not depend on cluster mass, whereas the dynamical 
lifetime of a cluster is strongly related to its mass. These differ- 
ences cause different schemes in the evolution of the luminos- 



4.1. Evolution of the mass function of star clusters 

In Fig. [6] we show mass spectra of three cluster subsamples 
differing just by the upper limit of ages of the clusters. The 
first sample contains the 49 youngest clusters with ages less 
than t y ~ 8 Myr (logf < 6.9). This is the lowest age limit 
where we have been able to construct a significant mass func- 
tion. Hereafter, we call this subset the youngest sample. The 
second subsample includes additional moderately young clus- 
ters and contains 207 entities (log f < 7.9). Finally, we consider 
all 440 clusters (log f < 9.5). To improve the stability of the so- 
lutions, we made re-binning on the mass scale requiring at least 
five clusters per bin, but keeping the step of AlogM f = 0.15 
whenever possible. 

All three distributions demonstrate similar general patterns 
that can be separated into two parts. The more massive part at 
log M c > log M* « 3.3 is almost linear on the logarithmic scales 
log rj, and log M c although the slope of the relation depends on 
the upper age limit. At lower masses, the CPDMF changes into a 
non-linear pattern. After reaching the maximum, the mass func- 
tion of each sample starts to decrease. For the youngest sample, 
the low mass portion of the CMF is rather flat and the turnover is 
less significant. For older clusters, the CMF maximum is shifted 
to lower masses. Choosing other subsamples of clusters by vary- 
ing the upper age limit from log t = 6.9 to log t — 9.5, we found 
that, at the age limit log t = 7.4, the CMF reaches a maximum at 
log M c ~ 2. If we add older and older clusters, this location does 
not change anymore, although the value of the CMF-maximum 
increases. 
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Table 1. Parameters of the CMF for the youngest sample. 



Fit 


Segment 


log 


,M C 


X 


log rf 






range 






S 


low mass 


1.70.. 


.3.30" 


-0.18 ±0.14 


-0.69 ±0.38 




high mass 


3.30" 


...4.90 


0.66 ±0.14 


2.16 ±0.56 




low mass 


1.70.. 


.3.37 2 » 


-0.18 ±0.14 




C 










-0.07 ± 0.60 




high mass 


3.37 2) 


...4.93 


0.66 ±0.14 





" logM r * = 3.30 

2) logM* from the solution of the continuity equation 



On the abscissa of Fig. [6] we observe different mass ranges 
for the different subsamples. The lower mass limit depends sig- 
nificantly on the cluster age. The younger the clusters, the higher 
is the lower mass limit (M""" ss 50M o at logf < 6.9 and 
M""" =a 5M at log t < 9.5). We consider this as a manifesta- 
tion of mass loss driven by the evolution of star clusters. On the 
other hand, the position of the bin of the highest masses seems to 
be independent of the sample considered. It varies from sample 
to sample only randomly, simply due to low number statistics of 
high-mass clusters. 

In Fig. [7] we show the evolution of the slope of the linear 
portion of the CMF. Again, one can split its evolutionary path 
into two parts. The younger clusters show a steadily increasing 
slope from^ x 0.65 at the upper limit t « 8 Myr to^- * 1.15 at 
t « 600 Myr. At older ages, the CMF slope practically does not 
change. 

The mass functions r], in Fig. [6] and their slopes x m Fig- 
describe cumulative mass distributions observed for clusters 
younger than t. Therefore, an increase in x with increasing t is a 
consequence of adding older clusters, the mass spectra of which 
are similar or slightly flatter than the CIMF. At the limit, the 
slope of r\, corresponds to the slope of the CPDMF representing 
a mixture of clusters of different ages. 

4.2. Initial mass function of star clusters 

In the youngest sample, even the least massive clusters seem to 
have not had enough time to evolve dynamically. Therefore, we 
assume in the following that the youngest sample gives a good 
approximation of the CIMF, and the youngest point of the evolu- 
tion curve in Fig.|7]is close to the slope of the CIMF. According 
to Fig. |6l we assume that Galactic clusters do form in a mass 
range from Mf" = 50 M to M n c mx = 3 • 10 5 M . 

The two segments of their CMF (logM e < logM* « 3.3, 
and logM e > logM*) can be scaled by a power law function 
with Eq. (HJ. The corresponding parameters are given in TableQ] 
(Fit S) where the fits are carried out separately for each of the 
two CMF segments. Again, since the CMF is based on clusters 
within the completeness limit, we consider the flattening J],{M C ) 
at log M c < 3.3 to be a real feature and not a signature of incom- 
pleteness. 

To convert Eq. ([TJ in a convenient form suitable for the com- 
plete range of cluster masses, we adjust the solutions making use 
of the CMF continuity condition, so the mass M* can be deter- 
mined where for both segments rj t (M*) = rj*. The results of these 
fits are given in Table Q] (Fit C). Obviously, the transformations 
do not affect the CMF slopes x- Applying these definitions, we 
obtain 



dN 



'7/ 



d log M c 



-x 



(2) 



For the construction of the CIMF we introduce a birth func- 
tion f3(M c , f) of star clusters as 



3 2 N 
dM c dt 



= p{M c ,t). 



Provided that the mass distribution of clusters at birth (CIMF) 
does not depend on time, we can factorize to get 

p(M c ,t) = aM c )v(t), 

where £(M C ) is the CIMF, and v(t) is the cluster formation rate 
(CFR). Provided tha t at t > tp the violen t events related to infant 
mortality (see e.g. Kroupaetal. 2001) have generally ceased 
and a cluster starts its post-natal life, one can assume that within 
a short time interval following this moment, i.e. between t o = 4 
Myr and t y - 8 Myr, the mass loss from a cluster due to two - 
body relaxation seems to be insignificant (lLamers et al.l l2005). 
and the relation between CMF and CIMF is straightforward: 



log e r], 



v(t)At M c ' 



where v(t) is the average cluster formation rate during the time 
At. Using Eq. (|2]i and the normalization condition 

J" &M c )dM c = \, 

M ,m„ 

we derive 

k iMrY" 



M* \M* 



(3) 



where or = 1 + and 

_ log e r]* 
v(t)At 

is the normalization factor. With the parameters from Table Q] 
one gets 



(4) 



0.82 + 0.14, 
1.66 + 0.14, 



for M"™ 
for M* 



< M c < M* 

< M c < M™ 



and M* - 2.3 ■ 10 3 M . The corresponding normalization factor 
is /c = 0.24. 

The two-segment CIMF can be understood in terms of the 
infant mortality process. The more massive clusters being able 
to survive the infant mortality and to keep the memory of then- 
parental conditions in the embedded phase show a power law 
distribution scaling with their nata l distribution. T his interpreta- 
tion is supported by the results of lLada et alj d 199 lb . who studied 
dense cores within the molecular cloud L 1630 and found that 
the mass distribution of dense cores can be fitted by a power low 
with an exponent of -1.6. The result is in remarkable agreement 
with -a = -1.66 we obtained for the CIMF at M c > 2.3 • 10 3 M Q . 

The low-mass segment of the CIMF probably consists for 
the major part of remnants of clusters damaged during the in- 
fant mortality phase and therefore lost memory of their recent 
past. Their mass distribution in the logarithmic scale resembles 
rather white noise. On the other hand, among young low-mass 
clusters, there may be a fraction that already had low masses at 
birth. Although the number of remnants and originally low-mass 
clusters is relatively large, their proportion is considerably less 
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Fig. 8. Evolution of the luminosity function 4>(Im v ) of Galactic 
open clusters. Different symbols mark samples with different up- 
per limits of cluster age. The filled circles are for clusters with 
logf < 6.9, stars show CLF for logf < 7.9, and crosses indi- 
cate the CPDLF based on all 440 clusters (logf < 9.5). The 
bars are Poisson errors. The arrow indicates the observed limit 
of integrated absolute magnitudes reached for open clusters in 
the LMC. The dotted vertical line at Im v = -8 separates the 
luminous part wher e the CLF defi ciency was observed for extra- 
galactic clusters bv lLarsenl d2002l) . 



than the one they had if they would follow the power law of the 
more massive clusters. 

Using the conclusions of a recent study of the influence of 
the gas expulsion phas e on the shape of the resulting CIMF 
(Par mentier et al.l 120 08). we infer that the empirical CIMF de- 
rived in the present paper is compatible with the scenario of a 
low, randomly distributed star formation efficiency with an aver- 
age of 20% and a dispersion of 3%. 

4.3. Evolution of luminosity function of star clusters 

The evolution of the luminosity function of Galactic star clus- 
ters is shown in Fig. [8] Similar to the mass function, the CLF 
apparently includes two regimes (although less clearly distin- 
guishable from each other than in the CMF case). The dominant 
feature is the linear regime of the CLF, occupying the bulk of 
brighter magnitudes up to Im v ~ -2.5, and a minor detail is a tail 
at fainter magnitudes. However, since the evolution of the lumi- 
nosity function is governed mainly by stellar evolution, and only 
at the end of cluster life by dynamical effects, there is no direct 
correspondence between absolute magnitude and mass of clus- 
ters. In most of the CLF bins of absolute magnitude, we find a 
mixture of masses and evolutionary states. Only at the extremes 
of the magnitude scale, there are clusters of more or less "pure" 
evolutionary status. As a result, the initial luminosity function of 
open clusters can be taken from the distribution of the youngest 
clusters (logf < 6.9) at Im v < -5.5 where the clusters belong 
to the massive portion of the CIMF (logM e > 3.3), and the tail 



of the CPDLF at Im v > -2 consists mainly of decaying clusters 
with masses less than * 1OOM i.e., after the CPDMF maximum. 

Contrary to the CMF, the luminosity function at -8.0 < 
Im v < -5.5 does not show any evolution of the slope with time. 
It seems that, with increasing age, the CLF is just moving along 
the ordinate due to the arrival of newborn clusters without chang- 
ing the shape and extends to fainter magnitudes due to the aging 
effect (fading of the brightest stars). 

5. Discussion of the results 

Based on this study on local Galactic clusters, we find that, in the 
immediate vicinity of the Sun, young (log f <6.9) open clusters 
are present within a wide range of masses. The mass range ex- 
tends from the highest masses - observed also in other galaxies 
- to a few tens of solar mass. This is a strong hint that the popu- 
lation in the solar neighbourhood is typical of spiral galaxies. In 
quite a natural way, cluster-like associations (see. Sect. |2| form 
the high-mass end of the CPDMF; in other words, they represent 
the same population of the Galactic disc as young open clusters. 
A comparison of our luminosity and mass functions of star clus- 
ters with those in other galaxies, which we discussed in SectfJ) 
strongly indicates a universal law of cluster formation, which is 
only weakly dependent on the environment. The issue of the uni- 
versality of the cluster initial mass function in diff erent galaxies 

is discussed in more detail bv lde Griis et al.l d2003[). 

Studying the CPDLF of extragalactic clusters, Larsenl d2002l) 
notes a depletion of the CPDLF at Im v < -8 in three galaxies. 
The author explains this effect by a poss ible increase in the steep- 
ness of the CPDLF. Giel es et all (|2006) confirm the existence of 
this precipitous drop in the cluster luminosity functions of the 
spiral galaxies NGC 6946 and M5 1 and explain this by assum- 
ing two-segment-, or Schechter CIMFs, in these galaxies. We 
observe the same feature in the luminosity function of Galactic 
clusters, too. Unfortunately, the ASCC-2.5 sample is restricted 
to the solar neighbourhood and does not include enough of the 
brightest clusters to be compared with the brightest limit of the 
extragalactic CPDLFs. As noted at the end of the last section, 
the parallel slope of the CPDLF for samples of different ages in 
the range -8 < Im v < -5 is a consequence of the continuity 
equation. A depletion at Im v < -8 suggests that no open clusters 
are born (or leave their parental cloud) with an absolute inte- 
grated magnitude much brighter than Im v = — 8. A caveat is the 
poor statistics (only a few extragalactic clusters and no Galactic 
ones). 

From the physical point of view, the observed deficiency 
might corresp ond to a gap caused b y so-called "type II embed- 
ded clusters" (Rroupa & Boilv 2002), which are massive enough 
to produce plenty of O stars providing fast gas expulsion and 
parent cluster decay. The low-mass limit of this cluster family 
sho uld be in the v i cinity of log M c = 4.5 in this case. 

Lada & Ladal d2003l) have compiled a catalogue of about 
100 embedded clusters within 2.4 kpc from the Sun. The sam- 
ple contains some optical objects and is partly ov erlapping with 
our da ta. Using models of the luminosity function. lLada & Ladal 
(2003) scaled the IR counts within the areas studied, estimated 
cluster masses, and constructed a mass function of embedded 
clusters. Typically, the clusters are distributed over a mass range 
from 50 to 1000 M and thus, related the low-mass part of our 
CIMF. Based on "reasonably complete subsamples", the mass 
distribution function M c dN/dlogM c increases approximately 
linear with M c . This indicates that the CMF of embedded clus- 
ters dN/dlogM c is nearly constant and agrees with low-mass 
CIMF (see Table[Qi. 
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The newly determined CIMF and CPDMF have major im- 
plications for the evolution of the Galactic disc population. The 
average mass of forming clusters can be calculated from Eq. ©: 



M, 



CIMF 



M"< 



I 



M""" 



M c ((M c )dM c 



With the parameters determined in Sect l4.2l we find that 

M c = 4.5 ■ 10 3 M G . The accuracy of this estimate depends on 
the accuracy of the CIMF parameters and, in particular, on the 
adopted integration limits. For example, by excluding the high- 
est mass bin (i.e., a shift of M" mx from 3 ■ 10 5 M to to 1O 4 M ) 



from consideration, one obtains M,, 



= 1.4- 1O J M 



The average mass of the presently observed clusters M t 
can be found from the CPDMF by its integration as 

log M'™* 

J M c ri(M c )d\ogM c 

CPDMF log Mf" 



CPDMF 



M. 



log M™ x 

/ r](M c )dlogM c 

log A/""" 

Taking the mass limits and the CPDMF from Fig. [6] we get 

CPDMF 

M c si 700 M Q . The accuracy of this value is estimated to 

CPDMF 

be about 30%. The value of M c coincides perfectly with 
what is usually assumed for open clusters. The average age of 
the presently observed clusters computed in the same way is 
equal to about 260 Myr. We stress that this estimate is based on 

... , . „ CIMF CPDMF 

purely empirical grounds. A comparison of M c and M c 
indicates that, during the first 260 Myr of its evolution, a typi- 
cal open cluster loses about 60%-80% of its initial mass. We 
consider this as evidence of the importance of this phase in the 
evolution of open clusters. The average mass-loss rate based on 
these numbers then ranges from 3 to 14 M Q Myr -1 . Since this 
estimate follows from observed masses of clusters, it includes 
both, mass loss due to stellar evolution and dynamical ejection 
of cluster stars. 

The average surface density of open clusters in the solar 
neighbourhood n c can be computed as 



log M™ 



Wc = jr 



tj(M c ) d log M c « 86 kpc 



This is about 25% low er than was estimate d in our study of 
the cluster population ( Pisk unov et al.l 120061) on the basis of a 
volume-limited sample of local clusters closer than 850 pc. This 
disagreement can be explained by our omitting of clusters with 
Iv > 8 in the present estimate. When taken into account, the 
estimates differ only by about 10%, which may result from the 
different methods of the sample extraction, i.e. the construction 
of volume-limited versus magnitude-limited samples. 

Let us estimate the input of open clusters into the Galactic 
disc population and consider the total mass of optical star clus- 
ters that were formed during the history of the Galactic disc. The 
surface density of the Galactic disc stars having passed through 
an open cluster phase can be written as 

CIMF 

Z = vTM c 

where T — 13 Gyr is the age of the Galactic disc, and v = 
0.4 kpc~ 2 Myr~' is the CFR determined via Eq. (0|. This value 



is determined as an average over the past 4 Myr and has to be 
compared with other determinations. For optical clusters, we 
find similar results in the literature: e.g. 0.45 kpc~ 2 Myr~' in 
Batti nelli & Capuzzo-Dolcettal dl991l) . and 0.25 kpc~ 2 Myr _1 in 
Elmegreen & Clemens! |l985). For embedded clusters, a much 



higher formation r ate of 2 - 4 kpc 2 Myr 1 is obtained by 
lLada&Ladal (12003). This may be regarded as evidence of a con- 
siderable infant mortality among embedded clusters. 

Assuming that the recent average CFR has not changed dur- 

CIMF 

ing the history of the Galactic disc one finds that for M c - 
4.5 ■ 1O 3 M 

I = 22M pc~ 2 . 

This amount corresponds to about 40% of the total surface den- 
sity of the Galactic disc in the sola r neighbourhood that, ac- 
cording to Holmberg & Flynn (2004) is 5 6 + 6 M Q pc~ 2 . From 
data on embedded clusters. lLada & Ladal d2003l) find X = 10 - 
3OM pc~ 2 . The value of S = 22M pc~ 2 obtained here is con- 
siderably higher than the previous estimates for the input of 
open clusters to the observed stellar popul ation of the Galactic 
disc that is qu oted as about 10% (see iMiller & Scald 1978; 
iPiskunov etall 120061) or even less than 10% (Wielen 1971). 
These rather low estimates can be explained by several situa- 
tions: they are based either on insufficient cluster statistics or on 
cluster masses from underestimated stellar counts, and on uncer- 
tain virial masses of open clusters or on arbitrarily chosen esti- 
mates of typical cluster contents. The newly derived value of E is 
in agree ment with indi rect evidence of binary population statis- 
tics (see iRroupall 1 995l) . This value significantly readdresses the 
problem of the main sources of the Galactic disc population and 
indicates open clusters to be one of the major factories of field 
stars, in line with the fact that a major fraction of star forma- 
tion seems to occur in embedded clusters (IZinnecker et all 1993; 
lLada & Lada1l200HlMegeath et al.ll2005l) . 



6. Summary and conclusions 

Using the integrated magnitudes computed from individual pho- 
tometry of cluster members available in the ASCC-2.5, we found 
that the sample of 648 clusters drawn from this catalogue is com- 
plete for apparent integrated magnitudes brighter than 8 mag, 
with 440 clusters above this completeness limit. The correspond- 
ing completeness distances depending on the brightness of clus- 
ters are typically between 1 and 3 kpc, although some prominent 
clusters can be seen up to distances of 10 kpc and more. 

The CPDLF constructed on the basis of the complete sample 
is observed in the range of Im v = [-10, -0.5] mag, i.e. about 
5 magnitudes deeper than observable in nearby galaxies. It in- 
creases from the brightest limit to a turnover at about Im v ~ -2.5 
and slowly falls down towards fainter magnitudes. We consider 
the turnover to be a real feature reflecting the behaviour of 
Galactic open clusters. At magnitudes brighter than the turnover 
position, the CPDLF is linear with a slope a = 0.41 + 0.01 (or 
a = 2.02+0.02), which is in perfect agreement with the observed 
slopes of cluster luminosity functions in nearby spiral galaxies, 
although these show a higher star formation activity. 

The masses of Galactic open clusters cover a range from sev- 
eral tens to hundreds of thousands of solar masses. The CPDMF 
mimics the general behaviour of the CPDLF: it shows a linear 
part at high masses (logM c > 2.5), a broad maximum between 
\ogM c = 1.5 and 2.5, and a slow decline towards lower masses. 
In spite of the present-day observations only allowing for a com- 
parison of the CPDMF data of extragalactic open clusters drawn 
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from model mass-luminosity relations, the agreement with pub- 
lished values is impressive. We find that for \ogM c > 2.5 the 
CPDMF can by fitted with high confidence by a power law with 
a slope x - 1-01 ± 0.04. We regard the observed coincidence 
of cluster mass function parameters derived for the Galactic and 
extragalactic clusters with different methods as indirect evidence 
that both methods of estimating cluster masses are calibrated 
consistently and give correct and reliable results. 

Inspection of cluster samples with different upper limits for 
age indicates that cluster mass spectra change with time. At ev- 
ery age the cluster mass function keeps the basic features of 
the CPDMF, i.e., a quasi-linear high-mass portion, and a non- 
linear portion at lower masses. With time, the slope of the lin- 
ear portion increases from x = 0.66 + 0.14 at logf < 6.9 to 
X = 1.13+0.08 atlogf < 8.5. Then the steepening of the CPDMF 
slows down and arrives at^ = 1.17 + 0.07 for logf < 9.5. The 
low-mass portion changes from an approximately flat distribu- 
tion at log t = 6.9 to a clearly non-linear behaviour displaying a 
broad maximum with a peak at about 100 M . 

We construct the CIMF from the data of a subsample of 
the 49 youngest clusters with logf < 6.9. The CIMF has a 
segmented structure with two power-law segments: the power- 
law slopes are a = 1.66 + 0.14 for logM r = 3.37 . . .4.93 and 
a = 0.82 + 0.14 for log M c = 1.7... 3.37. 

The luminosity function of open clusters does not show the 
same systematic steepening with age as does the mass function. 
It evolves via parallel-shifting along the ordinate without chang- 
ing the slope. Only the faintest segment of the CLF seems to 
depend on age and spreads towards the fainter magnitudes with 
age. The weak dependence on age and the presence of a mixture 
of clusters with different masses in every magnitude bin makes 
it difficult to determine the CIMF from the CLF. 

When comparing the average mass of the newly formed, 

youngest clusters M c - 4.5 • 10 3 M B with the average cluster 

— CPDMF 

mass from the whole sample (M c =2 700 M ), one finds, 
as expected, that clusters are typically much more massive at 
birth than later. Since this observation concerns "normal" optical 
clusters, one should acknowledge the importance of a continu- 
ous mass loss occurring in open clusters during their evolution 
with a rate around 3-14 M Myr~'. 

If this effect is not taken into account, one underestimates the 
role that open clusters play in supplying the Galactic disc with 
the products of their evolution: stars and stellar remnants, sub- 
stellar objects, gas and dust. Provided that the cluster formation 
history has not changed dramatically in the solar neighbourhood 
during the evolution of the Galactic disc, at least 40% of the 
observed surface density of the disc comes from open clusters. 
The other half is left for violent events (strong winds, HII zones, 
or supernova explosions) provoking the fast mass lo ss during 
the "in fant" phase o f cluster evolut i on pr oposed first by Tutukov 
dl978l) and late r by kroupa et alj d200l[). an d re cently consid- 
ered in detail by Baumgardt & Kroupa (2007j) and lWeidner et alj 
(2007b. 
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